Primordial non-Gaussianities of gravitational waves in the most general single-field 

inflation model 
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We completely clarify the feature of primordial non-Gaussianities of tensor perturbations in gener- 
alized G- inflation, i.e., the most general single- field inflation model with second order field equations. 
It is shown that the most general cubic action for the tensor perturbation (gravitational wave) hij 
is composed only of two contributions, one with two spacial derivatives and the other with one time 
derivative on each hij . The former is essentially identical to the cubic term that appears in Einstein 
gravity and predicts a squeezed shape, while the latter newly appears in the presence of the kinetic 
coupling to the Einstein tensor and predicts an equilateral shape. Thus, only two shapes appear 
in the graviton bispectrum of the most general single-field inflation model, which could open a new 
clue to the identification of inflationary gravitational waves in observations of cosmic microwave 
background anisotropies as well as direct gravitational wave detection experiments. 



Inflation, an accelerated expansion of the early Uni- 
verse caused by a scalar field called inflaton, is a quite 
promising paradigm of cosmology, and primordial pertur- 
bations generated from inflation are crucial clues to the 
yet unidentified inflationary model. From the proper- 
ties of the primordial perturbations such as the power 
spectra and spectral indices, we can extract informa- 
tion about the theory governing the inflaton dynamics. 
Among them the non-Gaussian signature in the cosmic 
microwave background (CMB) has been paid much at- 
tention in recent years, along with the great progress in 
precise cosmological observations. So far most of the lit- 
erature has focused upon non-Gaussianities of the scalar 
perturbations [l[ , as they are most directly connected to 
the CMB observations. Tensor perturbations 0, how- 
ever, are also generated during inflation, whose direct 
detection would be the most obvious evidence for infla- 
tion. When we try to detect tensor perturbations with 
the CMB measurements and/or with the direct detec- 
tion experiments, it is essentially important to remove 
the background (contamination) sources. For example, 
the B-mode polarizations are dominated by the lensing 
effects on relatively small scales. Astrophysical sources 
like white dwarf binaries could dominate the power spec- 
trum for a wide range of frequencies of the background 



gravitational waves. Thus, non-Gaussianities will be a 
key feature of the tensor perturbations [U, 0] as well as 
the scalar perturbations because they can help us to dis- 
criminate the inflationary signals from other contamina- 
tion sources even if the latter dominates the power spec- 
trum. For this purpose, we need to completely clarify 
the features of the non-Gaussianities of primordial tensor 
perturbations produced during inflation, which enables 
us to make templates for non-Gaussianities of primordial 
gravitational waves. 

In this Letter, we, for the first time, investigate 
the non-Gaussianities of primordial tensor perturbations 
based on the most general single-field inflation model, 
i.e., generalized G-inflation make a complete identi- 
fication of the shapes of bispectra, and explore the pos- 
sibility of large non-Gaussianities from the tensor sector. 

The Lagrangian for generalized G-inflation is the most 
general one that is composed of the metric g^ v , the 
scalar field <f>, and their arbitrary derivatives, and has 
the second-order field equations. The Lagrangian was 
first derived by Horndeski in 1974 [f|, and very recently 
it was rediscovered in a modern form as the general- 
ized Galileon |7jL i.e., the most general extension of 
the Galileon [aQ, in four dimensions. The generalized 
Galileon is described by the sum of the following four: 
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C 2 = K{4>, X), C 3 = -G 3 (<l>, X)n<t>, 

C 4 = G^,X)R + G 4X [(D0) 2 - (V, t V^) 2 ] , 

-3D0(V Al V^) 2 + 2(V A1 V^) 3 ], (1) 
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where K and GVs arc arbitrary functions of (f> and 
X := —{dcf)) 2 /2. Here we used the notation G iX for 
dGi/dX. The generalized Galilcon can be used as a 
framework to study the most general single-field inflation 
model. Generalized G-inflation contains novel models, 
as well as previously known models of single-field infla- 
tion such as standard canonical inflation, k- inflation [Io| , 
extended inflation |Tl|, R 2 inflation fl2|. new Higgs in- 
flation [l3| . and (minimal) G-inflation jl4j. The above 
four Lagrangians can even reproduce the non-minimal 
coupling to the Gauss-Bonnet term Q. 

In [5[, the background equations for generalized G- 
inflation is presented, and the most general quadratic ac- 
tions for tensor and scalar perturbations are determined, 
giving the power spectra of the primordial perturbations. 
The most general cubic action for scalar perturbations is 
worked out in flM Il6j . The curvature perturbation in 
generalized G-inflation is shown to be conserved on large 
scales at non- linear order in [l7| ■ We are going to present 
the most general cubic action for tensor perturbations to 
determine the possible tensor bispectrum arising from 
single-field inflation. 

The perturbed metric around a cosmological back- 
ground may be written as 



5oo 



1, 30, - 0, 9ij = a 2 (t) (e h ) i 



where 



( eh )ij = + hij + -jhikhkj + -^hikhkihij + 



(2) 



(3) 



and hij is a transverse and traceless tensor perturbation, 



dihi 







with repeated spatial indices summa- 



rized by dij. The perturbed metric defined in this way 
is convenient for calculating the action because we have 



We plug the metric ^ into the action 



5 = E 



-g£i 



(4) 



and expand it in terms of hij to get the quadratic and 
cubic actions. Only the Lagrangians that involve the 
curvature tensors and V^V„</>, i.e., £4 and £5, contribute 
to the quadratic and cubic actions. 

The quadratic action was already derived in Q : 



- / did 



3 xa 3 



Q T h 2 



where an ovcrdot denotes d/dt and 



F T ■= 2 G 4 -X(<j>G 



'5X 



G, 



5c6 



Qt 



G4 — IXG4X — X I HtfiG^x ~ G 



(5) 

(G) 
(7) 



From the action we see that the propagation speed of the 
gravitational waves is given by c/, := y 'Ft/Gt, which 
may differ from unity. In order for the system to be 
stable Ft > and Qt > are required. 



The linear perturbation equation derived from the ac- 
tion ([5]) can be solved in the Fourier space, 



f d 3 k ~ 
Mt ' x)= J (2^^*' k ) e * 



(8) 



To proceed further, it is convenient to introduce a new 
time coordinate defined by dy = c^dt/a. We assume 
for simplicity that the inflationary Universe may be ap- 
proximated by de Sitter spacetime and Ft and Qt are 
approximately constant. Using the normalized mode so- 
lution, 




yH^U-ky), 



(9) 



where H^ 2 is the Hankcl function, the quantized tensor 
perturbation is written as 

^•=$^[Vfce£ ) (k)o a (k)+# 



k e* W (-k)at(-k) ((10) 



(s) 

where is the polarization tensor with the helicity 



states s = ±2, satisfying e-^(k) = = kje\j'(k). Here 

we adopt the normalization such that e|^(k)e^- ^(k) = 
5 SS >. Choosing the phase of the polarization tensors ap- 



(k) 



propriately, we have the relations e*y (k) = eL s ^ 

efj (— k). The commutation relation for the creation 

and annihilation operators is given by [a s (k), er,(k')] = 
(27r) 3 £ ss /<5(k— k'). The 2-point function can now be com- 
puted as 

(/^■(k)Mk')) - (2 7 r) 3 ( 5( 3 )(k + k')^ J , fei (k), (11) 



>ij,kl 

where we introduced 
ny,fc{(k) = 



|Vfc| 2 IIij,fci(k), 



^eW(k)e^(k). 



(12) 



(13) 



The power spectrum Vh = (k /2tt )V%j,ij is given by 
1 H 2 



V h (k) 



2tt 2 Ftch 



(14) 



wherc y = —1/k corresponds to the time of the sound 
horizon exit . 

Having thus obtained the quadratic action and the so- 
lution to the linearized equation governing the inflation- 
ary gravitational waves, we now move on to the cubic 
action. The most general cubic action for tensor pertur- 
bations in the single-field context is obtained as 



did 3 x a 3 



Ft 

"4a 2 



X<j)G 5X 1 ; ; 
-^2 ^ij^jk^ki 



hikhji — -hijhki 



dkdihij 



,(15) 
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which is composed only of two contributions. Clearly, the 
term with one time derivative on each hij appears only 
if G$x 7^ 0. This term is absent in the case of Einstein 
gravity, non-minimal coupling to gravity (£3 = f(<p)R), 
and even in the case of new Higgs inflation which involves 
a non-standard kinetic term of the form G^ v d^ducj) [HI- 
However, in the presence of non-minimal coupling to the 
Gauss-Bonnet term, this term does not vanish [5[. The 
terms of the form h 2 d 2 h, where d represents a spatial 
derivative, is already present in the case of Einstein grav- 
ity, and even in the most general case only the overall nor- 
malization is generalized from the Planck mass squared 
Mpj to the function Tt- 

The 3-point function can be computed by employing 
the in-in formalism, 

(h nn (ki)^ M - a (k 2 )h i3h (k 3 )) 

= -i / dt'([/i njl (i,ki)/i 42j2 (t,k2)/ii3 J 3(t,k3),i7 int (t')]>, 



where to is some early time when the perturbation is well 
inside the sound horizon, t is a time several e-foldings 



after the sound horizon exit, and the interaction Hamil- 
tonian is 



flint (*) = - 



d 3 xa 3 



X<j>G 5X 1 ; ; 

12 "'jkibki 



(16) 



It will be convenient to introduce the non-Gaussian am- 
plitude Ai^izfeizjz defined by 

(h njl {ki)h i2 j 2 (k 2 )ft fs j 3 (k 3 )) 
= (2^^ + k 2 + ka)Pg ■ 

K l K 2 K 3 



We write A^^i 

now) 



A 



(now) 



A £S 2 hi 3 h> where 



A^ n ™K , o and A-4^} a i „ represent the contributions 

»Ul»2j2»3j3 2U122J223J3 r 

from the h 3 term and the h 2 d 2 h terms, respectively. Just 
for simplicity, here again the inflationary Universe is ap- 
proximated by de Sitter spacetime, which allows us to 
to compute the non-Gaussian amplitude assuming that 
X4>G§x — const and Tt — const. Each contribution is 
then found to be 



A 



(new) 



HXj>G 5X kiklki 

r^3 ,/"H K l ) L H 2 j 2 ,mn\K-2 )lh 3 j 3 ,nl ( K 3ji 



AGt K 3 



A n3\lj2^3]3 = A \ n *lil^( k l) n *2j2ji(k2) 



(18) 

5 perms of 1,2, 3 I , (19) 



where K = k\ + k 2 + k 3 and 



A(ki,k 2 ,k 3 ) := - — 



1 



1 k 3 ^ k * kj 



k\k 2 k 3 



(20) 



We see that the second contribution, which is present in 
the case of Einstein gravity, is independent of any func- 
tions in the Lagrangian, and hence for all the models of 
single-field inflation A^^} 2 j 2i3 j 3 coincides with the one 
in general relativity. (For this reason we associate this 
piece of the amplitude with the superscript "GR.") The 

size of the first contribution, -A^"^ j 2 i 3 j 3 , is crucially de- 
pendent on how (j> couples to gravity along the inflation- 
ary trajectory. Only those two amplitudes are sufficient 
to characterize the tensor bispectrum in the most general 
single-field inflation model. 

We are now in position to discuss whether or not large 
non-Gaussianitics can be obtained from A 



(new) 



»U1»2J2I3J3 ' 

Since Gt 3 HXipG^x, the ratio HXipG^x /Gt cannot be 
large in models with large HX<j)G$x- The only possibil- 
ity is that various terms in Gt are arranged to cancel each 
other to give Gt ~ 0. (Note, however, that Gt must be 
finite and positive.) This then yields J~s ~ —J~t, where 
J~s is the coefficient of (d() 2 in the quadratic action of the 



curvature perturbation £ and must be positive to avoid 
gradient instabilities (see Ref. [H). Therefore, models 
with small Gt tend to be unstable against either scalar 
or tensor perturbations. For this reason, generally speak- 
ing, it is rather non-trivial to get large non-Gaussianitics 
from the h 3 term, though one cannot completely deny 
the possibility of making both J-5 and Tt positive with 
the help of the functional degrees of freedom of our La- 
grangian. 

Let us turn to the two polarization modes, 



£<«>(k) :=^(k)e*«(k), 



(21) 



and consider their amplitudes j[ SlS2S3 of the bispectra 
(£ Sl (ki)£ S2 (k 2 )£ S3 (k 3 ))- The amplitude may be defined 
in an analogous way to Eq. (TlTl) . so that -4( n e W S ) 3 (GR) = 



tain 



(ki)ea ) (k 2 )e^(k 3 )^S- We thus ob 



* ( s 



(ncw),(GR) 



JS1S2S3 
(new) 

ASlS2S 3 

(GR) 



HX0G 5X k\k 2 k 2 

- 1 , 



4Gt 



K 3 



sis 2 s 3 
(now) 



A(k 1 ,k 2 ,h)F^(k 1 



^ (fci,fc 2 ,fc 3 ),(22) 
k 2 ,k 3 ), (23) 
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AUUUi/k, Jfcs/*i)/(*a/*iX*3/*i) 



k 2 /k 




as a func- 



FIG. 1: ^Jio(l,*a/Jfei ) *3/*i)(*3/*i) _1 (*a/*i)~ 1 » 
tion of kijk\ and kzjk\. The plot is normalized to unity for 
equilateral configurations k%jk\ = fe/fei = 1. 




FIG. 2: J; 2 /Ai, feAOtfe/fci)" 1 ^/*!)" 1 as a func- 



tion of k^/k\ and k^/ki. The plot is normalized to unity for 
equilateral configurations fo/fci = fe/fei = 1. 



where we defined 



A' 3 



A 3 -4^fc, 2 fc i -4*1*2*3 



. (24) 



*(+++(*!, fc 2) * 3 ) := A 2 A+++(fc 1 ,fc 2 ,fc 3 )/2, and 



^&), ( GR)(*i.*a.*3) := f M,(GR)( fc i' fc 2,-fe). Since 
our theory accommodates no parity violation, we have 

F — und P H — F+H — 

r (ncw),(GR) ~~ r (new),(GR) r (new),(GR) ~~ r (new),(GR)' 

The non-Gaussian amplitudes and A^t^ are 

plotted in Figs. [T] and [2] One sees that the amplitude of 
the new contribution peaks in the equilateral configura- 
tion, while the "GR" contribution becomes largest in the 
squeezed limit. This gives a clear distinction between 
the two different contributions, and the two character- 
istic shapes would be helpful to discriminate the infla- 
tionary gravitational waves from those produced by other 
sources. The other correlation functions such as the — ++ 
one are subdominant relative to the + + + one because 



for equilateral configurations i^cw) 



F +- + 
(GR) 



(GR.) 



(new) (new) 



F+++/81, and 



(GR) . 



in the squeezed limit, k 3 

U2 



-k'i/2 and 



p-++ 

(GR) 



0, one has A^p| 
^3 ■ 



F, 



(GR) 



In this Letter, we have clarified primordial non- 
Gaussianities of tensor perturbations arising from the 
most general single-field inflation model with second- 
order field equations, and have found that they are com- 
pletely determined by two different contributions: -A( new ) 
and ^4(gr)- Our results provide at least two distinctive 
features to test the framework of generalized G-inflation 
based on the graviton non-Gaussianities. Firstly, *4( ne w) 
is a unique feature of the kinetic coupling term G5 . Any 
detection of this type of bispectrum, no matter large 
or small, would unambiguously indicate the existence of 
non- vanishing G§x, at least in the Galilcon framework. 
Secondly, the contribution -4(gr) is a fixed and univer- 
sal feature for single-field inflation models which are all 
within the generalized G-inflation framework. It is im- 
possible to enhance/suppress this contribution in gener- 
alized G-inflation models. In other words, any detec- 
tion of enhancement / suppression of this contribution to 
the graviton bispectrum would imply new physics be- 
yond generalized G-inflation and/or other astrophysical 
sources. The two contributions are clearly distinguish- 
able according to their shapes of non-Gaussianities. 
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